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THE MOVEMENT OF WATER DROPLETS IN CLOUDS 
AROUND THE NOSE OF AN ATMOSPHERIC RESEARCH AIRCRAFT 

P. Feuillebois, M. Scibilia 
Polyphasic Flow Team 
CNRS* Aerothermal Laboratory 


INTRODUCTION 

The underlying technical problem presented by this study 
is the movement of water droplets in clouds around the nose 
of an aircraft. In practice, this consists of collecting 
samples of suspended droplets which are representative of the 
cloud, for the purpose of future meteorological investigations. 
It is natural to think that larger drops will break over it. 

The proportion of the different droplet sizes, and their res- 
pective concentrations would therefore be different on each 
point around the aircraft nose. The position of the droplet 
sensor tube, in front of the aircraft nose, should therefore 
influence the samples collected. In this report, we will 
calculate the droplet concentrations, to determine more pre- 
cisely what they are. 

From the theoretical standpoint, we tried to evaluate the 
dynamic interactions between droplets and air flows. We will 
show later in this report that the interactions between drop- 
lets in a cloud are negligeable. In other words, each particle 
behaves (for the calculation of the forces exerted upon it) as 
if it were the only one in the air. Moreover, we will show 
that, in most cases, the flow of air carrying droplets is not 
influenced on the average by the presence of these droplets. 
These two results then make it possible to study separately 
the trajectory of each droplet after having calculated the 
flow of air assumed to be dry. To calculate the concentrations, 
we will introduce the Langragian method. The concentration is 
calculated directly according to two close trajectories. 

♦National Center for Scientific Research. 
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Experiments are performed at the National Astronomy 
and Geophysics Institute [I.N.A.G.] on an atmospheric re- 
search aircraft. The technical data were reported to us by 
GRIBKOFF, in charge of devising these experiments. 

The first position of the sensor tube is schematized in 
figure 1. The numerical data are the following: 

-aircraft speed: v m = 60 m/s, 

-pressure: p = 800 to 1100 mb, 

-temperature: T^ = -10°C to +30°C /2 

-saturating steam concentration: 10 to 15 g/m" 5 of air 
-water droplets: 

•diameter: d = 1 to 50 p 

3 

.concentration: m = 1 g/m of air on the average 

3 

3 to 4 g/m of air maximum 
.number of droplets: n = 200 to 500 per cm^ of air. 

1 - CALCULATION OF DIMENSIONLESS NUMBERS AND DISCUSSION OF THE /3 
PHYSICAL CONDITIONS 

1.1. REYNOLDS NUMBER 

To determine the rate of flow about the nose of an air- 
craft (assumed to be hemispherical, of radius 0.7 m) , let us 
calculate the Reynolds number. 

To accomplish this, the air viscosity must be evaluated. 

Dry Air : According to SCHLICHTING [1] we have: 
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Influence of Humidity: 


According to [2], the saturating water steam pressure is: 
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The molar concentration is maximum for 30°C and therefore: 

P - 2 

X = _s = 4.3x.O 

P 


Such a small molar concentration does not influence the 
viscosity [ 2 ] . 


Consistency of data pertaining to humidity : 


/4 


At 10°C , we have a number of air moles per m of air: 


n 


air 22.4 


WOO x 273 = 43 


283 


This gives a water weight per m of air: 


m H 0 = 18 n n = 18 x n . = 18 x 1.252x10 2 x 43 = 9.1 g 

2 H2 u 3 

which is consistent with the datum 10 to 15 g of saturating 
steam per m^ of air. 


Influence of Water Droplets: 


According to the data, m = 1 g of water/m of air on the 
average. The concentration of droplets per unit of volume is 

g 3 

therefore, with p = 10 g/m Cwater density) : 

ir 
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m n „-6 

c = 10 5 6 g/m 3 10 

The viscosity of a diluted suspension is given in a 
first approximation by Einstein' s formula, on the assumption 
that the water droplets follow the air flow: 


y = y . 
air 


(1 + J). 


Since c = 10 6 , we see that the variation of y is negli- 
geable. in regard to the density of the mixture, the concen- 
tration of droplets is about: 


in la 

— = jtjqq - 1%, and is therefore negligeable also. 

dir 


Finally, the water droplets do not influence the kinem- 
atic viscosity v = — . 


In conclusion, the kinematic viscosity of humid air load- 
ed with water droplets in clouds is virtually equal to that of 
dry air. 


The Reynolds number based on the speed, v^ of the aircraft 
and on the radius R, of the hemispheric nose of the aircraft 
has the following values: 
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According to SCHLICHTING, the flow about a sphere placed 
in a flow with little turbulence becomes turbulent for: 

5 

Re is greater than or equal to 1.92x10 
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Moreover, the atmosphere, where the turbulence is charac- /5 
terized by large vortices, may be considered on the sphere 
scale as having little turbulence. 

The turbulent flow about a sphere is characterized by 
the boundary layer separation situated downstream from the 
midship bulkhead. The flow is turbulent in this region. 

However, downstream, the flow is virtually that of an inviscid 
flow. 

This is the case here where 

Re is greater than or equal to 2.6x10^ 

The turbulent region thrown downstream has no influence 
on the flow about the hemisphere, which is an inviscid flow. 

1.2. POTENTIAL FLOW ABOUT THE SPHERE 


The inviscid fluid flow about the sphere is a potential 
flow, which may be expressed in cylindrical coordinates (r,z) 
(figure 1).. 

v z 3 

CO R 

^ = V oo z ” 2 ~~2 2 3/2 

* (r + z ) 7 

The components of the relative air speed relative to 
the aircraft are, according to r,z, respectively: 
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For the sake of convenience, we will employ dimension- 
less magnitudes using R as the reference length, and as 
the reference speed. Let us put: 


c “ i 


v - JL . 3 _CP 
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At the measuring point shown in figure 1 (tube end) , the 
reduced cylindrical coordinates are: 


/ 6 
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On this point, the reduced speed components calculated 
are equal to: 

|V | = 0.2626 
P M 

|V P | = -0.6817 
^ M 

The angle which the speed forms with the z (or £) axis is: 
8 = Arctg & | = 21° 


V 


K 


and the reduced speed modulus is: 


V ■ // 2 + v _ 2 ■ 0 , 73 OS 
P 5 • 
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We see that the measuring point is within a region 
where the speed is fairly inclined over the axis, and where 
the modulus is appreciably smaller than the value at infinity. 

As nothing leads us to suppose that the droplets follow 
the fluid flow, the movement of the droplets should be stu- 
died in more detail. 

1.3. STOKES NUMBER 

This number characterizes the movement of particles with 
respect to the movement of air. It is expressed: 



-t is the deceleration time of a droplet. On the assump- 
tion that the droplet is subjected to the Stokes drag, this 
time is about: 

4 

. _ i 1 E (at the radius of the droplet) : 

‘p 6 it > p 

. _i 

-t is the characteristic time of the air flow about the 
nose of an aircraft: 


We therefore calculate: 

SK = | Re ^ (f ) 2 
P 

With the values calculated for Re, and the values for 
air density (assuming an inviscid gas) : 
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Consistency of the Data With the Droplet Sizes Under Consideration 

Above, we calculated that the droplet density concentration 

3 

is on the average c = 10 for m = 1 g water/m . At the maxim- 
al value m = 3 to 4 g water/m 3 of air therefore c = 3 to 4x10 
Another datum is the number of droplets 

n = 200 to 500/cm 3 


On the basis of c and n we may deduce the droplet diameter 
(already given; the data are superabundant) : 
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We calculate the values of d: 
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the droplet diameter is within the interval 15-30 yim (con- 
sistent with the data d = 1 to 50 yra) . The most common in- 
terval found is 15-30 ym. 


Given the preceding, let us calculate the typical val- 
ues ofthe Stokes number: 
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Physically, a low Stokes number corresponds to droplets 
which have a tendency to easily follow the air movement. This 
is the case of droplets whose size is smaller than or equal to 
10 ym. A high Stokes number corresponds to droplets which con- 
tinue a straight forward movement, and are not at all influenced 
by the movement of air. For the majority of the droplets con- 
sidered here, a Stokes number of order 1 corresponds to an 



intermediary state in which the droplets follow the fluid 
flow "slightly". This is theoretically the most complicated 
case where it is necessary to calculate the trajectory of 
each type of droplet individually. 

2 - EQUATIONS FOR THE SUSPENSION OF OF DROPLETS IN ITS MOVE- - /9 

MENT AND NUMERICAL CALCULATIONS 

2.1. DRAG FORCE ON THE DROPLETS 


In the expression of the Stokes number, we have assumed 
that one droplet is subjected to the Stokes drag force. This 
formula is valid for droplets that are fairly small such that 
the Reynolds number of the relative flow about each droplet 
remains low. We will see below that this condition remains 
verified, except for the largest particles of 50 ym for which 
we shall give a correction for the drag. 

A second condition for applying the Stokes formula is 
that the droplets behave dynamically as solid spheres, i.e. 
that the movement of liquid in the droplet may be disregarded. 
According to LEVICH [4] , this is verified for diameters between 
0.5 ym and 10 urn. Without more date, we will still keep this 
assumption for larger diameters. 

The Stokes force is valid for an isolated particle. How- 
ever for a suspended particle, the interactions between par- 
ticles may be involved. However the correction to be made is 
of the order of c, and since c - 10 , it is discardable here. 

The fact that c is small implies that the droplets be- 
have as specific forces (their volume is negligeable) . 

2.2. EQUATIONS GOVERNING SUSPENSION 


With the preceding assumptions, we will express for the 
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mixture of fluid flow and droplets a system of movement 
equations, according to MARBLE. These equations are gen- 
eralized here to N types of droplets; they are expressed 
as dimensionless (see Feuillebois' thesis): 

-Fluid continuity: £ v~0 

-Fluid momentum: 

N ($ - V) 

(V . $)V ■ - + — . V 2 V + Z f, - P i '~ . 

e _ i«i Sk i 


-Continuity of type i particles: V*. 


(f. V .) 
1 pi 


0 


-Momentum of type i particles: 



^)V - . = (V - V . )/SK, 

pi pi 3 


with the following dimensionless variable: V fluid velocity, 

P fluid pressure, and , for type i particles, V ., velocity, 

pi 

f^ concentration per unit of mass (mass of particles per 
unit of fluid mass), Sk^ Stokes number. 

It is interesting to observe that the term f ^ (V^-V/SK^ 

for £he effect of the particles on the flow is negligeable. 

Actually, if we calculate the concentration per unit of mass 

f. from the numerical data, we find: f. = m./p where m. = 1 
l 3 i 3 1 K 1 

to 4 g/m of air, p = 0.92 to 1.45 kg/m . Therefore, 

_3 

f^ is less than or equal to 4 . 3 x 10 

The terra f . (V . - \b/Sk. could be nondiscardable if SK. 

1 pi 1 1 

were small at the same time as f., i.e. if the small particles 
were the most abundant in the mass. However, we have seen 
earlier that the most abundant particles are in the range 
15-30 ym. 

The interaction term f^ (.V^-X^/Sk^ is therefore always 
^discardable . 
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Consequently/ the air flow about the aircraft nose is 
not modified by the presence of particles. The system of 
equations cross-checks and it is simply necessary to solve 
independently the system pertaining to each type of particle. 
The equation for the momentum of the particles is also ex- 
pressed in Lagrangain form, and in cylindrical coordinates: 


dV 


dT 


J2_ = 


V_ - V 


PP 


sk 


dV 

dT 


PS 


V_ - V 


Sk 


£i 


d/dT is the derivative with respect to time calculated on 
the basis of one particle in its movement. 


T is the reduced time 


T = 


t v 


00 


R 


Furthermore, and to simplify the annotations, we have 

eliminated index i relative to each type of particle. To 

calculate the droplet trajectory, the values given above of 

V ,V for the inviscid fluid flow and the expressions of 
P S 

the particle velocity: 


V = dp_ 
P P dT 




If we combine these equations, we obtain the following 
system to solve: 



dT 


( 3 n ^ - r —) — 

<n 2 - C 2 ) V2 7 Sic 


12 - v 
dt pp 




i » 15, s v 

dt PC 


The initial conditions to be used express that at upstream 
infinity the droplets have the velocity of air: 

T -*• - ® s p * p„ i C -* ro 

v ♦ o i v >1 : 

P P P C 

For the numerical calculation, we replaced this condi- 
tion at infinity by the following condition: 


T = 0: p = p Q ; C -> 


V -> 0; V - -1; 
r P 


which proved to give the required accuracy. 

2.3. CALCULATIONS OF THE TRAJECTORIES 

The integration of the differential system was achieved 
using the DREBS sub-program available at CIRCE. . The sub-pro- 
gram, in FORTRAN, performes dual-precision calculations. 

The droplets .trajectories were calculated, then sketched 
for four cases: 


Sk =0.01; 0.1; 0.5; 1 


/12 


13 



characteristic of the droplet sizes studied. 


In figures 2 to 5, we see that the smallest droplets 
(Sk = 0.01 by-pass the aircraft nose, whereas the largest 
(Sk = 0.5 Sk = 1) will hit against the aircraft nose, as 
could be expected, according to the definition of the Stokes 
number (see above) . 

It is also remarkable to note that these numerical re- 
sults coincide with the theoretical calculation of MICHAEL 
[5] who predicts that the particles collide with the sphere 
for Sk is greater than 0.12. 

The trajectories were calculated in detail near the 
opening of the measuring probe, or near: 

P M = 0.571 ± 0.011 
C M * 1 ' 18 

2.4 CALCULATION OP THE PARTICLES CONCENTRATION 


We replace the equation for the particles continuity 
in a form more adaptable to the present numerical calcula- 
tion. 

For each type of particle (fixed Sk) , let us consider 
two very close trajectories, each differing from the other 
by the initial value of p (figure 6) . 

P _ P 

P = P D + Ap o 

We express that the total particles flux across the 
axisymetric surface made up of two surfaces of particles 
streams, and of end crowns, is zero: 



We have : 


f V . AS + f ^ . AS. = 0 

P ° 


Using a little geometry (figure 6), gives us: 


AS = 2irp Ah 


Ah = Sp cosa 


<5p = Ap + As sina 


where : 


therefore: 
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PC 
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The conservation of flux is therefore expressed: 




V Ap 
_E_ P ° 


p (V Ap + v AC)V 
P . C 


where : f 


V - / V 2 + V 2 

p p o p ; 
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We can therefore calculate on the basis of two close 
projecteries, the ratio of the concentration at the meas- 
uring point f to the concentration at infinity f Q . 

The error with such a method may be considerable if the 
trajectories are not evaluated with sufficient accuracy/ as 
the study below will show. 

2.5. CALCULATING THE ERROR ON THE VALUE OF CONCENTRATION F 


f - — - 


V P 


p (V Ap + V A;)v 


We will note 6X the error range on the quantity X. The 
relative error on F is calculated: 


if - jgo , 60 , 

F 0. p V p Ap o 


«v • 6[V Ap + V Ac] 


p o 


VpAP + V p AC 
P 


p Q is the initial value set for p, therefore 5p Q = 0; 
P Q + Ap is another intial value, therefore <5(Ap)= 0 


/ v 1 + v 2 
p p P <. 

|v I + |V IdV 






_Q_ 


A(V p Ap + V p Ac] - V p 6(Ap) + Ap «V p + V p | 5 ( Ac ) + | 4; | 6V 
P . C p p 

Let us take the example of calculations made for Sk = 1, 
the tolerance being 10 5 . We obtained the following results: 
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Checking the calculation accuracy on a computer: 
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If we consider the values : 


p^ - 0.53 > T - 8.9 » p --0,56510 > C * 1.2310 ; V - 0.074882 ; V - 0,86767 
° Pp Pc 
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Ad - 0,01 j A<5“ 0,01058 » AC ■ -0,00080. 
K o 


we obtain: 


SE 

p 


1,49.10 


The most important contribution comes from the error 

in As , due to the error in [V Ap + V AsJ 

P Pp 

2.6. VALUES CALCULATED FOR CONCENTRATION 
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The results are given in figures 7 and 8. 
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2.7. CORRECTIONS ON THE DRAG FOR LARGE PARTICLES 


When the Reynolds number of the relative flow with re- 
spect to one particle is low, we may use Stoke's force. 

This is no longer true if the Reyonlds number, based on the 
diameter, becomes greater than 1. (see the figure extracted 
SOO's book [6] . 



The relative Reynolds number based on the diameter 

2a | v - v I 
Re = P 

e d 


v 


may be evaluated according to the results already calculated 
using S'toke ' s force on each droplet. The number Re^ will in- 
crease with the size of the particles, for which the error 
will greater. 

Let us take the example: 2a = 50 ym (for these droplets, 
Sk - 1).. 
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According to the calculation results, the velocities 
at the measuring point are: 

V - 0,2626 ; V - -0,6817 

P C 

V ■ 0,U7489 } V ■ -0,86765 

p p A 


Therefore 


and 


|v p - v| - / (v - v p )« + tv - v ? )“ 


2a v «> l V p ~ V i 50X10"* X 50 [ Vp - V| 


, Re v- V.. 


1.5x10 


-s 


52.84 


For this Reynolds number, the coefficient of the /17 

sphere's drag {see "Standard drag curve") is, according to 
the figure, about 2.45 times that of Stokes law. There is 
no good empirical formula giving the coefficient of drag for 
intermediary Reynolds numbers. The least undesirable one still 
seems to be that of KLYACHKO Isee FUCHS, 1966 [7]) and PUTMAN 
(1961), cited by RUDINGER 18]. 


c d ■ (1 + I Re a /3) 

24 

where Stoke 's coefficient, — is multiplied by an empirical 

K 

e d 

factor. For R =52.84, this factor is equal to 3.35. The 
e d 


curve borrowed by SOO gave the factor 2.45, but we should 
rather have more confidence in the formula above. RUDINGER 
actually indicated, through several comparisons with exper- 
ience, that its accuracy is within a few %. 


To convert the formulas, we simply have to multiply Stoke ’s 

1 2/3 

coefficient 6 fray by (1 + — Re^ ) which is the same as divid- 
ing the Stokes number 


4/3 TT 


SK = 


6iray 


-P by (1 + | 


Re a 2/3 ) 
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In this factor Re^ is variable: 


_ 2a v 

Re a = |v„ - v 

V 1 o 1 


We introduce the new dimensionless numbers according 

to Reynolds number Re = RV^/v of the aircraft and the adi- 

2a 

mensional size of droplet D = — . Therefore: 


Re, = D Re V - V 
d P 


The system to be integrated is therefore expressed: 


cv 


" - v ) (1+7 [D Re (V n - V) ] ) 

P 


V: , JL 


dT ''P 'p/ “ 6 '' J ' Sk 


dV 


DT 


^ - ( y c - V p } (1 + j [D Re (V p - V)] 2/j ) — 


Sk 


with: 


V ■ - 
P 2 


-e_s_ 


(p 2 + z 2 ) h 


/’B 


V . -i - 1 ItL 1 - 2g» ) 
5 7p‘ + r 2 ) 


V » / v * + v 4 
p c 


v/"v*v 

p . p c 


and 


dT V p 


^ - v 
dT V 
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D and R values: 

e 

The diameter 2a of the droplets varies from 1 to 60 ym. 
Therefore: 

D = ^ + ppy to - 6 - °Q 1 y- 6 = 1.43xl0~ 6 to 8.571xl0~ 5 

Re varies from 3.48xl0 6 (at -10°C) to 2.6xl0 6 (at 30°C) . 
Therefore 3.76 < D Re < 298 

Values Selected For the Calculation : 

We will select the mean values corresponding to: 

T = 0°C; p = 1013 mb 

Therefore 

Re 


Sk 

by setting D 
Therefore : 

-2 —1 

We select for Sk the values 10 ; 10 ; 0.5; 1. 

At the temperature and pressure selected, these values 
correspond to the droplets diameters: 


/19 

] 
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3 CONCLUSION FOR THE AIRCRAFT WITH HEMISPHERIC NOSE 


/20 


Even with the drag correction, the values of the con- 
centration F remain close to one. This means that, to the 
nearest 3%, the values f measured in the tube will actually 

be the values f in the cloud, 
o 

The 3% deviation is the maximum, corresponding to drop- 
lets of about 40 ym. For smaller or larger droplets, this 
deviation is less. 
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- DROPLET TRAJECTORIES FOR Sk = 0.1 

- DROPLET TRAJECTORIES FOR Sk = 0.5 

- DROPLET TRAJECTORIES FOR Sk = 1 
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Fig 1. Diagram of the aircraft nose. 
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Droplet trajectories 
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Fig. 5. Droplet trjaectories for Sk 





Fig. 7. Concentration values as a function of the droplets 
diameter . 


Relative increase of concentration per unit of mass at the 
measuring point with respect to the concentration per unit 
of mass at infinity as a function of droplet diameter in 
normal conditions. 

(.The segments represent error ranges in the numerical calculation) 
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